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Abstract. The temperature dependence of static dielectric susceptibility of a system 
with strongly coupled fluctuating dipoles is calculated within a self consistent mean 
fluctuation field approximation. Results arc qualitatively in good agreement with, a 
quantum paraelectric, SrTiC>3 in the low temperature regime. We identify this system 
as a gaped quantum paraelectric and suggest a possible experimental realization of a 
quantum critical paraelectric through the application of hydrostatic pressure or doping 
by impurity. 
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1. Introduction and summary 

Phase transitions in displacive systems can not be described by Ising type Hamiltonian, 
which is usually invoked for a system going through order-disorder transition. Behavior 
of these systems are essentially governed by collective oscillations of coupled dipoles 
and the phase transition is described by softening of the corresponding optical mode 
due to thermal fluctuation [TJ [2]. There are materials like SrTi03, KTaC>3 which 
are supposed to show a displacive transition similar to that occurs BaTi03, but fail 
to do so. In fact there is softening of optical phonon modes in these materials also, 
but that does not lead to a transition even at very low temperature [3j H]. Instead, 
there is a strong enhancement of low temperature dielectric constant. For example, 
in SrTiOa the static dielectric susceptibility has a very high saturated value ((9(10 4 )) 
at low temperature (~ 10K) followed by a Curie regime (10 to 100 K) and a long tail 
thereafter. It seems that the classical soft mode concept is insufficient to describe various 
aspects of low temperature behavior. Explanation of the unusual behavior of dielectric 
susceptibility and the stability in the low temperature paraelectric phase has been a 
long standing puzzle. Long ago Barrett [5] gave a semi-phenomenological theory, which 
essentially recasts the Curie- Weiss formula with a replacement of temperature T there, 
by average energy, thereby the inverse of dielectric susceptibility could be written as, 
X -1 oc Ti coth(Ti/T) — T c , where T c is classically calculated critical temperature and Ti 
is a quantum scale (9(/?/mass). This theory, in the high temperature limit, reproduces 
the Curie Weiss law. To match experimental data in SrTiOs the Barretts' formula has 
been found inadequate as one single constant quantum scale 7\ can not trace the full 
curve. The formula has since been modified in various ways, for example, by introducing 
an extra exponent [6], that is, by writing x~ l a s (T\ coth(Ti/T) — T c )~ u , and by making 
T\ temperature dependent with an extra scale [7] , to take care of various "anomalies" , 
for example the one near 40K. There has been a proposal of attributing this extra energy 
scale to the structural transition which occurs at 110K [8j. These proposals either follow 
an order parameter expansion similar to the Landau expansion or some modifications 
thereof, hence they do not introduce any new microscopic description. In SrTiOa and 
KTaOs there is no ordering, therefore the low temperature regime where the dielectric 
constant is enhanced, the physics is dominated by fluctuations of relevant microscopic 
degrees of freedom rather than their averages. In this article we analyze the fluctuation 
in such systems within a self consistent mean fluctuation field approximation. There are 
mainly two parameters, namely the anharmonicity parameter and the effective stiffness. 
The zero point or quantum fluctuation will be dominant when stiffness is small. The 
qualitative behavior of susceptibility is reproduced as well as a new insight gained into 
the quantum critical behavior of such systems. A mismatch in theory and experiment 
curves for the dielectric constant at high temperature can be attributed to the effect 
of structural transition which occurs at higher temperature (i.e. at 110 K in SrTiOs), 
such discrepancy is irrelevant for the following discussion which refers mainly to the low 
temperature regime. 
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2. Mean Field Analysis 

The low temperature physics of SrTiOs is dominated by fluctuations of Ti ions [9]. The 
Hamiltonian for such ions is modeled in terms of local quartic oscillators coupled with 
a nearest neighbor harmonic interaction [10] , 

-^vu lUv . (1) 

l K ' IV 

The constants A and v are assumed to be positive and mass taken as unity. This 
Hamiltonian describes two local minima with a nearest neighbor coupling v. For 
\v\ « \uq\ and Uq <C 0, it mimics a two state Ising system with Gaussian fluctuations 
around one of the local minima. When \v\ ~ \&q\, there is a possibility of large tunneling 
between these minima. In this regime the system has to be described in terms of its 
collective behavior. Such system is called displacive system and the limit \v\ — * \uq\ is 
called Displacive limit. In momentum space, 

H = Y1 \ p2 i + \ J2( u o- v5 cosqa)u q u_g+^X ^ u qi u q2 u q3 u_ qi _ q2 _ q:i (2) 

Here 5 is the coordination number and a is lattice spacing. Now with p q = u q = 
—luouq in the kinetic energy term, finally the Hamiltonian within the quasi harmonic 
approximation, i.e., 

u qi u -qi (3) 



where 



is given by, 



where 



- = £<r« ff (o)u_,(o + )) I (4 ) 



H = 2 ~ uj2 ) u i u ~<i ( 5 ) 

i 



uj q =ujQ— vS cos qa + 3\a ~ cUq — v + vSa 2 q 2 + 3\a (6) 

is the renormalized frequency for small q (such truncation is quite justified for low 
temperature properties of a near critical system.). We are interested in the paraelectric 
phase of the system, that is, where < u >= 0. Since the system is at low temperature 
and the dielectric constant has an enhanced value, < u 2 > need not vanish, however. The 
purpose of present work is to present a self consistent calculation of < u 2 > in classical 
as well as in the quantum regime. The susceptibility, which is related to < u 2 >, is 
essentially the phonon propagator, 

x (q,n) = — ^ -, u n = 2nnT. (7) 

{iuj n Y - ui 
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With ujq given by equation^ we have a self consistent equation, 



/3 



-T 

ft ^ 



q,n 



q,n 



The solution of this equation will give <j(T) which in its asymptotic forms reduces to, 
a = V — ~ / alq q 2 — (Classical) (9) 



alq q 2 1 = (Quantum). (10) 

a/Wq — v + vog z + 3Ao" 



To go into details of temperature dependence of susceptibility, we need to define 
some physically interesting dimensionless parameters as, A = — (u 2 . — v)/u)q, <j c = 
— {uJq — u)/3A, Tj — Hj (2uoa c ) (h is taken as unity in this article), as such A ~ a c ~ rj~ l . 
The parameter A describes the effective stiffness for collective modes at harmonic level. 
The strength of coupling between various modes near q = is determined by cr" 1 
while the parameter 77 tells us about the vicinity to the quantum limit in the system. 
Introducing normalized temperature x = Tjmuj\a c and using the previously defined 
parameters, we rearrange the equation ([6]) as follows 

4 = ^ + A (^-l) (11) 

and 

^ = y^coth^V (12) 

a c uj„ KujqxJ 



A self-consistent solution of these equations will give, 

X^O^oc A(--l). (13) 

For large enough A the system shows classical behavior, that is a ~ T from equation (Q. 
The mode coupling would give corrections higher order in temperature, and T c would be 
proportional to A. On the other hand as A become smaller and r\ becomes larger and the 
system move towards the quantum domain. When A or T c becomes identically zero we 
have quantum critical point. At this point the zero temperature static long wavelength 
susceptibility also diverges. Interestingly the A = or u> 2 = v limit is the displacive 
limit, well known in the structural transition literature. A non-self-consistent estimate, 
with a temperature dependent momentum cut-off (q ma x ~ T), tells that a starts from 
a constant value in the low temperature side and then follows a T 2 behavior in high 
temperature (up-to Debye temperature) side. Such a non-self-consistent estimate gives 
quite correct result when the system is far away from the quantum critical point i.e. 
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|(a;o — v)/3\o~\ << 1. At quantum critical point, an estimation of the self-consistent 
correction is also ~ T 2 . The Barret formula can not reproduce this result. That formula 
is essentially outcome of quantum fluctuations of in a single mode theory, which would 
fail near the quantum critical point as many modes and their coupling would dominate 
the behavior of system there. This necessitates a self-consistent calculation for quantum 
paraelectrics near its quantum critical point. From figure(l), we learn that the high value 
of static dielectric susceptibility of SrTiC>3 is attributed to the smallness of the parameter 
A(= 0.003). This motivates us to treat this system to be near quantum critical point. 
The static dielectric susceptibility data of SrTiC>3 remind us of the behavior of itinerant 
Fermionic systems near quantum phase transition point and fluctuation regime around 
that. There the (staggered) magnetic susceptibility diverges for (anti-)ferromagnetic 
transition as the coupling constant crosses a critical value [TT]. The case of SrTiOs is 
similar to that of liquid Helium-3 [12], where the magnetic susceptibility gets enhanced, 
as large as ten times. denendinEr imon nressure. from its free Fermionic value. 




Temperature (Kelvin) 



Figure 1. Numerical solution shows saturation in Static susceptibility (in units of 
10 ) vs Temperature curve. This curve is in good agreement with Muller's experiment 
in low temperature side, with = 1, A = 0.003, 77 = 1/A,q max = 0.1 and at the 
end x an d T are rescaled with 0.4/ A and 30A respectively.The lower curve is the non- 
selfconsistent fit with the same parameters as the upper one but with rescalling of \ 
and T by 9.5 and 100 respectively. 



3. Quantum Criticality and Hydrostatic Pressure at QCP 

After the above identification we now focus on theoretical aspects of quantum criticality 
in ferroelectric systems. At zero temperature zero mode fluctuation is the most 
dominant. From the equation ( fl2l) it is clear that zero temperature fluctuation a s 
is given by 

^0 V 
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0.5 

I Gapped Quantum Paraelectric 

Figure 2. Schematic phase diagram of a typical quantum paraelectric system 

Writing a s /a c as y we get, 

y 3 -y 2 -a = 0, a = ^ (15) 

This tells if a = then cr s = is a solution, which is the classical limit. For a ^ all 
solutions become non-zero and since y 3 — y 2 = a > that means y > 1. Thus cr s > (To, 
moreover the cr s increases as "a" increases. The meaning of quantum criticality, in terms 
of y, is y — > 1. In that limit the zero temperature properties will show a scaling behavior 
with (inverse of Correlation Length) x _1 ~ A(y — 1) ~ r/ 2 ~ (1 — v/uJq) 2 . If we define 
a quantum scale £q = (1 — v/uJq)~^, then x ~ £q 4 . The point ?/ = 1 is essentially the 
point where effective stiffness (A) changes sign. In the regime A > self consistency 
in fluctuation breaks down, system seeks ordering and hence an expansion about the 
non-zero < u > is required. In this case transition temperature T c ~ A2. On the other 
hand, in A < regime the system can not have any ordering and its behavior has to 
described by self consistent fluctuations. There is a characteristic temperature (cross 
over temperature in modern parlance [13]) T* ~ A^ which demarcates the boundary 
between the low temperature gapped quantum paraelectric behavior and the classical 
behavior(figure(2)). In case of SrTiC>3, the plateau in the susceptibility vs temperature 
curve is the signature of gapped quantum paraelectric behavior. There is no transition in 
this system. But there is a crossover from low temperature quantum to high temperature 
classical behavior at the crossover temperature T* (~ 10K ). This is exactly the 
temperature where plateau ends and the susceptibility curve stars following a Curie 
behavior (figure(l) and (3)). One can now hope to reach at A = through tuning some 
parameters like pressure, impurity etc. The width of this plateau regime vanishes at 
this point and the system becomes quantum critical. At this point thermodynamics will 
be described by power laws in temperature (e. g. x(0? O)" 1 ~ T~ 2 ) and the system will 
show some non trivial dynamics. The later is beyond the scope of the present work. It 
is quite evident here as the controlling factor v/moj 2 strongly depends upon structural 
aspects and hence this quantum-ness in SrTiOs can be properly understood through 
some intrinsic mechanism which give rise to such large tunneling . 

A good possibility of exploring the physics near such quantum critical point is 
through application of hydrostatic pressure. Such a technique is already used in case of 
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ferroelectrics and quantum paraelectrics long ago [2] and more recently [T3] in different 
contexts. We found that those experimental results can be discussed more interestingly 
as is done in the context of itinerant magnetic system recently[16]. Application of 
hydrostatic pressure will couple to optical mode via its coupling to the acoustic mode. 
In this case the starting Hamiltonian takes the form 

H = \ J d( l[p 2 g + M ~ v5 cosga) u q u- q ] + j J Uf =1 (dqiU qi )5(^2qi) 

+ g J dk dq e(k) u q u k _ q + y J dqe 2 (q) - pe(0) (16) 

Here last three terms are results of applications of pressure, in lowest possible order. 
The parameter il g n couples strain fields to unit cell displacement related to optic mode 
and U K" is the force constant for harmonic acoustic phonons, and the last term shows 
the coupling of the hydrostatic pressure "p" to the static strain with some unit strength. 
Now if the pressure is strong enough e has a minima at e = e(0) and is given by 

e(Q)=p/K. (17) 

Integrating out strain field, we get an effective Hamiltonian 



H= dq 



\pq + 2 M + 9P ~ VS COS 1 a ) U Q U ~q 



1 



with renormalized coupling constant of quartic term 

*«=(*-ir) (19) 

Again we write a self consistent equation for paraelectric fluctuations as, 

a = j d d q^ coth (^) (20) 



Where 

uj 2 (q) = 3AA(1 + p/po) + v5q 2 a 2 /2 + 3X R a and p = 3KAX (21) 

9 

Up to this point result is just renormalization of the factor A as A(l + p/po) and it 
becomes an experimentally controllable parameter. And the behavior of susceptibility 
at different values of A is shown in the figure. In this proposal we assume the positivity 
of \r. Otherwise transition will be first order and the scaling behavior will not be 
valid. In real life one can try to induce the effect of negative pressure required in these 
systems to achieve QCP through some homogeneous effects of non-polar impurity. But 
in either case nature of the transition can be modified because of strain coupling or 
disorder respectively. Analysis of such transition in these materials will be discussed in 
an upcoming paper [TTj. 
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Temperature 



Figure 3. Temperature variation of Susceptibility at different values of A and the 
log-log plot of the same 

4. Discussion 

We have shown that a mean fluctuation field theory within a quasi harmonic 
approximation reproduces the low temperature behavior of susceptibility of a quantum 
paraelectric quite well. The short range model studied here is justified since only 
transverse optical modes are involved in the ferro-electric fluctuations. In presence 
of a long range force longitudinal mode becomes stiff and only transverse modes can get 
soft. Presence of long range dipolar forces can induce a certain amount of anisotropy 
in the transverse modes, which can certainly change the critical behavior, however, 
only with a fairly large value of dipolar contribution to anisotropy in the quadratic 
term [18]. We are not aware of the first principle results for anisotropy parameters in 
case of SrTiOs or KTaC>3. However, the first principle calculations support our choice 
of parameter for the effective stiffness. Compared to BaTiOa it is about twenty times 
smaller (Table V in ref [19]) for SrTiOs, which makes it more near the quantum domain. 
On the other hand the lattice induced anisotropy in the quartic term is of the same order 
of magnitude and it would not play a key role in distinguishing the low temperature 
behavior in these systems. We leave discussions on anisotropy dependence for the future 
work and stick to isotropic short range model. It is also clear that there is no need to 
introduce "anomalous" regime as proposed earlier. That proposal might be due to the 
insistence on comparing experimental results with Barrett's formula and its extensions. 
The experimental behavior is well accounted for in the quantum region and at high 
temperature the susceptibility smoothly crosses over to the classical behavior. Here we 
have focused more on physics of low temperature behavior than the exact calculation of 
various properties. Thus the structural aspects and anisotropy effects are not attempted 
here. 
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